Discrepancy between sub-critical and fast rupture roughness: a cumulant analysis 
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We study the roughness of a crack interface in a sheet of paper. We distinguish between slow 
(sub-critical) and fast crack growth regimes. We show that the fracture roughness is different in 
the two regimes using a new method based on a multifractal formalism recently developed in the 
turbulence literature QJ. Deviations from monofractality also appear to be different in both regimes. 
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Since the early description of rough fractures as self- 
affine surfaces [2|, the existence of universal roughness 
exponents has been strongly debated 0- There are now 
many experimental evidences for a non-unique value of 
the roughness exponent of fracture surfaces. Different 
exponents can be found due to the anisotropy of the frac- 
turation process , the anisotropy of the material struc- 
ture Q or anomalous scaling related to finite-size effects 
0. A recent observation also shows that, in rupture of 
paper, the crack interface can not be described with a 
single roughness exponent and would be multifractal Q • 

Most of the time, roughness exponents appear inde- 
pendent of crack velocity. For rather slow velocities, 
no effect of the velocity on roughness has been ob- 
served in plexiglas (v = 10~ 7 — 5 10~ 5 m.s^ 1 ), glass (v — 
10" 9 - 5 10~ 8 m.s^ 1 ), intermetallic alloys (v = 10" 8 - 
510~ 5 m.s _1 ) or sandstone (v — 10~ 4 — 10 _2 m.s^ 1 ) 
0. On the contrary, in dynamic fracture of plexiglas 
(v > 600 m.s -1 ~ 0.45x Rayleigh wave speed), varia- 
tions of the roughness exponents with the velocity have 
been observed [tjj- As pointed out in a recent review |l0j| . 
there is a lack of experimental studies concerning the in- 
fluence of fracture kinetics on roughening. 

In this Letter, we study the roughness of a crack in- 
terface in a sheet of paper ^lj for which multifractality 
was observed Q. We show that it is crucial to take into 
account the non-stationarity of crack growth. Schemat- 
ically, crack growth occurs in two different regimes: a 
sub-critical regime where the growth is slow (v = 10~ 5 — 
10- 2 m.s _1 ) and a regime where the crack growth is fast 
(v ~ 300 m.s -1 ). We show that the fracture roughness 
is different in the two regimes. Various method of analy- 
sis of fracture roughness have been proposed Here, 
we introduce a new method based on a multifractal for- 
malism recently developed in the turbulence literature 
pj. This formalism allows us to measure reliably scal- 
ing exponents in fracture and to quantify very precisely 
deviations from monofractal behavior. 

Experiments. We recall briefly the experimental set-up 
described in . We use paper which is a bi-dimensional 
brittle material with a quasi-linear elastic stress-strain 
response until rupture. Samples are fax paper sheets of 
size 24 x 21 cm 2 manufactured by Alrey. Each sample 
has an initial crack at its center and is loaded in a ten- 



sile machine with a constant force F perpendicular to 
the crack direction (mode I). The stress intensity factor 
K(L) oc Fy/~L, where L is the crack length, determines 
the stress magnitude near the crack tip and is the control 
parameter of crack growth. For a given initial length Li, 
sub-critical crack growth is obtained by choosing F so 
that K(Lj) is smaller than a critical threshold K c corre- 
sponding to the material toughness. During an experi- 
ment, L increases, and so does K(L). It will make the 
crack accelerate until reaching the critical length L c for 
which K(L C ) = K c and above which a sudden transition 
to fast crack propagation occurs. Using a high speed and 
high resolution camera (Photron Ultima 1024), we have 
determined which part of the post mortem crack interface 
corresponds to slow or fast growth, and measure the ve- 
locity of the crack in each one. In the sub-critical regime, 
the velocity ranges from 10 -5 to 10 m.s . Recording 
at 4000 fps, we find a crack velocity about 300 m.s -1 in 
the fast regime. Note that there are four to seven decades 
between the two growth regime velocities. 





k s ( x ) 

| c 


b 


a 


b 


c „, r 













FIG. 1: Digitized post mortem sample and corresponding ex- 
tracted front showing three stages: (a) initial crack (Li = 
2 cm) , (b) sub-critical growth, (c) fast growth. 

Crack profiles. Post mortem samples are digitized with 
a commercial scanner at 1600 dpi, which corresponds to a 
pixel size oq — 16 /im close to the typical diameter of cel- 
lulose fibers. Clearly, below this scale, there is no inter- 
esting information about the scaling behaviour of cracks. 
On Fig. ^ we show an example of a digitized sample 
compared with the extracted crack front s(x). We distin- 
guish between different stages corresponding to: (a) the 
initial crack, (b) the sub-critical crack growth and (c) the 
fast crack growth. We have digitized 51 fractured sam- 
ples, obtained with different forces (200 N, 230 N, 250 N, 
280 N) and initial crack sizes (1cm, 2 cm). Since the ini- 
tial crack is centered, each sample give rise to two fronts. 
Thus, we have a total of 102 independent fronts of about 
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10 3 points for slow growth and 10 4 points for fast growth. 

Scale invariance. Let {s(x),x £ R} be a signal as a 
function of a coordinate x. Scale invariance of s means 
that there is no characteristic scales in the signal. The 
scaling properties of s can be characterized by introduc- 
ing a multiresolution coefficient T[s](a;, a) defined at scale 
a and position x. Scale invariance implies that the g-th 
order moments of the multiresolution coefficient follow a 
power law with exponents C,(q): 

M q = (\T[s](x,aW)^C\a\^ , (1) 

where the bracket denotes the average over the x space. 
The increments over a scale a, T[s](x, a) = s(x+a) — s(x), 
are standard multiresolution coefficients and the corre- 
sponding moments are the structure functions It 
has been shown that a more general framework for defin- 
ing multiresolution coefficients is the wavelet transform 
[iSEIil- When the signal s follows Eq. Q with ((q) 
proportional to q, the signal is monofractal. The com- 
plete analysis of the deviations of £(g) from monofractal- 
ity can be made through the multifractal formalism. 

Multifractal analysis. Multifractal formalism is based 
on the mathematical definition of a local singularity 
exponent h(x). Using the multiresolution coefficient 
T[s](x, a), we write at each position x [lij ]: 

T[ S ](x,a)~C\a\ h W , (2) 

where h{x) is the Holder exponent or local roughness 
exponent describing how singular the signal is at position 
x: the larger h(x), the smoother s(x). The statistical 
distribution of the Holder exponents is quantified by the 
singularity spectrum D(h) defined as |l4l llq: 

D(h) = d H {x\h(x) = h} . (3) 

where dn is the Haussdorf dimension. The probability to 
observe an exponent h at scale a is then proportional to 
fl i-_D(/i)_ Thus, the £(g) spectrum can be related to the 
singularity spectrum D(h) by a Legendre transform, i.e. 
C(q) — min^l + qh — D(h)). When s(x) is monofractal, 
h(x) is a constant H independent of x, D(H) — 1 and 
C(<z) — qH is proportional to q. Conversely, if s(x) is mul- 
tifractal, h{x) takes different values at different positions 
x, and £(g) is not proportional to q. 

In practice, the values ((q) are obtained by fitting- 
straight lines (when power law behavior is observed) on 
log-log plots of M q versus scale a for different moment or- 
ders q. The singularity spectrum D(h) is then deducted 
from To verify if £(g) differs from a linear monofrac- 
tal behavior, one needs to obtain £(g) for a large range 
of q values and then proceed to fit the ((q) curve (for 
instance, 1 < q < 8 in p|, see a ^ so 0D- 

While this has been the traditional way of estimating 
C{q), an alternate method introduced recently in the tur- 
bulence literature Q, 0] , involves only a few straight line 



fits (as low as 3) while still accurately estimating the non- 
linear behavior of the C(g) spectrum. To summarize this 
method, we start with the general expansion Q]: 

oo n 

\nM q = VCntojSL, (4) 

n=l 

where C„(a) are the cumulants of Q a = hi\T[s](x,a)\. 
One can demonstrate that the first three cumulants are 
the mean, standard deviation and skewness of Q a : 

Cl(o) = (Q a ), 

C 2 (a) = (Q 2 a )~(Q a ) 2 , (5) 
C 3 (a) = (Q 3 a )-3{Q 2 a ){Q a ) + {Q a ) 3 . 

Identifying the first derivative of Eq. and of the loga- 
rithm of Eq. with respect to ln(a), one finds: 

((q) = Cl q + c 2 q 2 /2\ + c 3 q 3 /3l + --- (6) 

where Ci = dCj(a)/dln(a) are constants in the case of 
scale invariant signals. It turns out that the average value 
of the Holder exponent is (h) = C\ and its variance (h 2 ) — 
(h) 2 = — C2/ln(a). When the multiresolution coefficient 
T[s](x, a) has Gaussian statistics (for example, when s is 
a Brownian motion), C*2(a) = 7r 2 /8 and C2 = 0. 

The above developments imply that ((q) can be 
estimated from linear regressions of the cumulants 
C n (a) vs ln(a) For a monofractal signal, c n = 0, Vn > 
2 and only one linear regression is needed. For a multi- 
fractal signal, a quadratic approximation requires 
only two linear regressions. In comparison with the stan- 
dard method based on the g-th order moments, the effi- 
ciency of the cumulant method becomes apparent. 

In the following, we will plot only results obtained us- 
ing the increments for T[s](x, a). Quantitative compari- 
son with other methods will be presented afterwards. 
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FIG. 2: (a) log 10 Mi and (b) C\j ln(10) versus scale calculated 
by the increments method for five different couples (F,Li) 
during fast crack growth. 

Influence of force and initial length. On Fig.0 we plot 
log 10 Mi and C\j ln(10) versus log 10 (a/ao) for five differ- 
ent couples of values (F, Li) during the fast growth stage. 
We see no dependence of Mi and C\ on the force or the 
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initial crack size. The same independence is observed for 
the slow growth stage. These observations are indepen- 
dent of the definition chosen for T[s](x,a) which means 
that our analysis is robust. In order to improve statis- 
tics, we will average the scaling laws obtained for data 
with different forces and initial crack lengths. Wc also 
notice that close to the discretization scale ao the slope 
becomes close to unity. This effect can be attributed to 
the discreteness of the signal In the following, we 

will concentrate on larger scales a > 4ao (— maximum 
fiber diameter) where this effect can be neglected. 
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FIG. 3: (a), (b): Mi and M2 for slow (circles) and fast 
(squares) crack growth. Statistical errors of the mean (= 
standard deviation /vN with TV = 102) are smaller than 
symbol size, (c), (d): Corresponding local slope with an error 
comparable to the symbol size. 

First and second order moments. Fig. EUa) and (b) 
show log- log plots of Mi and M2 versus a/ao for slow 
and fast growth. The curves are slightly different for slow 
and fast crack growth. Another important observation is 
that they are not perfectly straight lines. To magnify the 
difference, we plot the local slope of those curves (esti- 
mated as the mean slope on a half-decade window) on 
Fig. |3{c) and (d) . One can see an evolution of the slopes 
with scale, in a similar fashion for fast or slow cracks. 
However, we observe a systematic difference between the 
slope values which is almost independent of the scale but 
depends on the moment order. Also, assuming that the 
moments are approximately linear, it is not possible to 
conclude for a roughness exponent value since the typical 
slopes C(l) and C(2) are such that C(2) ^ 2((1). 

First and second order cumulants. On Fig. 0Ja) and 
(b), we plot Ci(a) and 02(a) divided by ln(10) versus 
log 10 (a/ao) for slow and fast crack growth. C\(a) is al- 
most linear, but has the same qualitative behavior than 
Mi or M2 when looking at the local slope (see Fig.Qfc)). 
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FIG. 4: (a), (b): Ci and C2 versus scale with statistical error 
bars. In (b), the horizontal solid line corresponds to a signal 
with Gaussian statistics, (c), (d): Corresponding local slope 
of cumulants with an error comparable to the symbol size. 



We will note the local slope cf (a) for the fast regime and 
cf (a) for the slow one. The separation between fast and 
slow crack growth is clearer on C\ than on Mi and M2, 
and the difference in local slopes Aci(a) = cf (a) — cf (a) 
is also more regular. The local extremum observed for 
a/a ~ 10 1 ' 5 (a ~ 500 ^m) in Fig. ^c) suggests the 
existence of a characteristic scale for both regimes. In 
contrast, 02(a) is hardly a linear function. We note that 
there is a range of scales (10 < a/ao < 10 2 ) for which 
the value of C2 is close to the theoretical value for a sig- 
nal with Gaussian statistics (see also [2l|). However, in 
this range, the values of the local slope can be non-zero, 
especially for the slow regime (see Fig. Eld)). 

Scaling laws. From the various plots, we can already 
conclude that it is not so easy to find a range of scales 
for which clear scaling laws are observed. For a/ao J$ 10, 
the slope of C\(a) is changing a lot because we start 
to feel the discretization effect previously discussed [ioj . 
We clearly see that for a/a > 10 2 , the slope of C\(a) 
is again changing significantly and seems to go towards 
0.5. This can be attributed to a change of statistics at 
large scales. Thus, if a scaling law exists, it is observed 
mainly at intermediate scales where a small oscillation 
of the local slope is observed. The same conclusion can 
be reached by looking at 02(a). At large scale, 02(a) 
decreases to values smaller than the Gaussian value that 
have no physical meaning. At small scales, C2 is very 
sensitive to discretization effects. At intermediate scales 
(10 < a/ao < 10 2 ), O2 appears quasi-constant. If one 
tries to estimate the slope C2, one finds that C2 is very 
close to zero for the fast crack growth and between —0.08 
and —0.02 for the slow part. 
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TABLE I: Scaling exponents ci (= (h)) for fast and slow crack 
growth, and their difference using various methods. 



Method 


Fast growth 


Slow growth 


Difference 


SF 


0.64 ± 0.02 


0.70 ± 0.02 


0.06 ± 0.01 


CWT 


0.65 ± 0.02 


0.73 ± 0.02 


0.07 ± 0.01 


WTMM 


0.64 ± 0.01 


0.70 ± 0.02 


0.06 ± 0.01 



length, is an open issue. 

We thank S. Ciliberto for fruitful discussions. This 
work was funded with grant ANR-05-JCJC-0121-01. 



We notice that the extremum value of c 2 (~ —0.02) 
observed at intermediate scales for slow crack growth is 
close to the one found in where both regimes were 
mixed. It is still difficult to confirm multiscaling since 
c 2 can be considered constant only over a very small 
range of scales. Nevertheless, the non-zero values ob- 
served for c 2 are important to consider since they will lead 
to C(2) 7^ 2^(1). Indeed, stopping at the second order, 
Eq. © predicts: C(l) = ci +c 2 /2, and C(2) = 2(c x + c 2 ). 
Thus, C(l) and £(2) are both influenced by the mean 
value of the Holder exponent and its variance. In order 
to give meaningful information about the scaling proper- 
ties of the crack, it is best to estimate directly ci from 
C\ since, in the case of scale invariant signals, c\ = (h). 

Scaling exponents. We have estimated the scaling ex- 
ponents in the range of scales for which the scaling laws 
are reasonably good. We took a/ao = 10 for the lower 
cut-off value, and for the upper cut-off value the point 
where C 2 cross the Gaussian value (a/eto = 10 2 for fast 
cracks, a/aa = 10 1 ' for slow cracks). In this scale range, 
we measure the mean and standard deviation of cf , cf 
and Aci . We show in Table [I] the mean values ± the 
standard deviation using various methods (SF: structure 
functions, CWT : Continuous Wavelet Transform 1J| 
and WTMM : Wavelet Transform Modulus Maxima [lg 
using I s * derivative of Gaussian). All methods give a dif- 
ference between the two growth regimes. For instance, 
structure functions give a difference of 0.06 ± 0.01 with a 
roughness exponent of 0.64±0.02 for the fast regime and 
0.70 ± 0.02 for the slow regime. Note that these values 
are consistent with the ones in the literature 0, El • The 
standard deviation is smaller on the difference than on 
the scaling exponents themselves due to the regular shift 
between local slopes seen on Fig. 0Jc) . 

Conclusion. Whatever are the exact scaling properties 
of cracks in paper, we find that the roughness changes 
when the rupture mechanism goes from sub-critical to 
fast crack growth. A decrease of the roughness exponent 
due to dynamical instabilities was previously observed in 
the case of fast crack growth 9] . It would be interesting 
to understand if the smaller roughness in the fast regime 
comes from a similar effect. Nevertheless, roughness in 
the slow regime is probably controlled more by the ma- 
terial disorder than by dynamical effects. There seems 
to be a characteristic scale a ~ 500 ^m in both growth 
regimes where the local slope of moments or cumulants 
reach an extremum value. Whether this scale is related 
to some internal structure of paper, for instance the fiber 
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